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Abstract 

Thousands of long intergenic non-coding RNAs (lincRNAs) are encoded by the mammalian genome. However, the function 
of most of these lincRNAs has not been identified in vivo. Here, we demonstrate a role for a novel lincRNA, linc-MYI-l, in adult 
fast-type myofiber specialization. Fast myosin heavy chain {MYl-i) genes and linc-IVIYI-i share a common enhancer, located in 
the fast MV/-/ gene locus and regulated by Sixl homeoproteins. //nc-MVH in nuclei of fast-type myofibers prevents slow-type 
and enhances fast-type gene expression. Functional fast-sarcomeric unit formation is achieved by the coordinate expression 
of fast /MVHs and iinc-MYI-i, under the control of a common Six-bound enhancer. 
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introduction 

Adult skeletal muscles are composed of slow and fast myofiber 
subtypes which selectively express the genes required for their 
specific contraction activity and metabolic properties [1-4]. These 
properties are acquired at the end of fetal development and during 
the neonatal period, when mixed skeletal myofibers expressing a 
panel of embryonic, fast and slow genes develop a specific slow or 
fast phenotype. The formation of efficient fast sarcomeric units, 
and Ca''"'' cycling and excitation/contraction/relaxation coupling 
in fast-myofibers, is achieved through the coordinate control of fast 
Mjihs and associated fast sarcomeric genes (including Tnnt3, Tnni2, 
Tnnc2, Atp2al and Pmlb) [2,4] . Myofibers can be classified by their 
MYH expression profile: slow-type myofibers in mice express 
MYH7 (also known as MYHCI, P or slow), and fast-myofibers 
express MYH2 (MYHCIIA), MYHl (MYHCIIX) or MYH4 
(MYHCIIB). Fast Mjh genes found in developmental and adult 
stages iMyh3, Myh2, Myhl, Myh4, Myh8 and Myhl3) are organized 
as a cluster within a 300 kb region on mouse chromosome 11 [5]. 
The spatio-temporal expression of one specific fast Myh gene at the 
Myh locus is reminiscent of the organization and expression of 
Globin genes at the beta-globiii locus [6]. However, we are yet to 
investigate potential enhancers or the Myh locus control region 
(LCR) that could be responsible for sequential and specific Myh 
gene expression in myofibers. The coordination of fast-type and 
slow-type gene expression in fast myofibers is not currently 
understood. Distinct intramyofibrillary calcium transients, evoked 
by slow tonic motor neuron firing, induce a cascade of 
downstream signaling involving Calcineurin and CamK. This 
results in the activation of selective transcription activators and 



repressors in slow myofibres. However, the signaling pathways 
operating in distinct Myh2, Myhl and Myh4 myofiber subtypes, 
which coordinate the activation of the other fast-type genes and 
the repression of slow-type genes, is less well understood [1]. Better 
knowledge of the mechanisms controlling muscle specialization 
and plasticity is important to enable the understanding and 
modulation of muscle adaptations in pathophysiological condi- 
tions. 

Six homeoproteins are major myogenic transcription factors 
which directly bind to DNA sequences (called MEF3s) to control 
myogenesis [7,8] and the genesis of fast-type myofibers during 
embryogenesis [9,10]. In adult skeletal muscle, Sixl accumulates 
at a higher level in the nuclei of adult fast myofibers than in of slow 
myofibers. Forced expression of Six 1 and its Eya 1 cofactor in slow 
myofibers causes adult slow-twitch oxidative fibers toward a fast- 
twitch glycolytic phenotype [11]. Animals with a Sixl KO present 
severe muscle hypoplasia and die at birth [12]. This prevents the in 
vivo analysis of the adult phenotype and the ability to investigate 
the direct or indirect involvement of Sixl in the spatio-temporal 
control of the expression of genes in the fast Myh cluster. 

The mammalian genome encodes thousands of long intergenic 
non-coding RNAs (lincRNAs) which have multiple functions 
[13,14]. Some accumulate in the cytoplam as miRNAs decoys 
[15,16]. Others accumulate in the nucleus and participate to gene 
regulation through chromatin remodeling and epigenetic modifi- 
cations [14,17,18]. Here, they may act as cis [19] or trans [20] 
transcriptional activators, as transcriptional repressors [21,22] or 
through DNA-RNA triplex formation [23,24]. 

In this study we identify a new lincRNA, linc-MYH, and the 
mechanism of its control of adult muscle fast fiber-type specification 
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Author Summary 

Adult skeletal muscles are classified into fast-type and 
slow-type, which display different resistance to muscle 
atrophy and metabolic protection against obesity. We 
identify in this manuscript a new mechanism controlling in 
vivo adult muscle fiber-type specification implicating a 
long intergenic non-coding RNA, linc-IVIYH. We demon- 
strate a three-element genetic partnership, where an 
enhancer under the control of the myogenic homeopro- 
tein Sixl functions as a regulatory hub to control fibre 
phenotype. In this partnership, the enhancer controls 
positively the expression of both the adjacent fast myosin 
heavy chain (IVIYH) gene cluster and of linc-MYH. linc-MYI-l 
is present only in adult fast type skeletal myofibers and 
controls their phenotype by suppressing slow-type gene 
expression. The regulation of iinc-IVIYI-i could provide a 
lead for new therapeutic approaches or drug develop- 
ment. 

in vivo. We demonstrate a three-element genetic partnership, 
where an enhancer element under the control of the myogenic 
homeoprotein Sixl functions as a regulatory hub to control fiber 
phenotype. In this partnership, the enhancer positively controls 
the expression of both the adjacent fast Myli gene cluster and 
linc-MYH, suppressing slow-type gene expression and facilitating 
fast fiber-type specialization. 

Results 

Sixl binds directly to a newly identified enhancer of the 
Myh genes cluster 

Our previous studies suggested that Sixl could be directly 
involved in the control of the expression of fast Myh genes, since 
higher levels of this transcription factor accumulate in the nuclei of 
adult fast myofibers than in slow myofibers [1 1] . To investigate 
how Sixl could control the expression of fast Myh isoforms, we 
used computational analysis to locate MEF3 sites at the fast Myh 
locus (see Materials and Methods). Six clustered MEF3 sites are 
conserved across human, rat and mouse genomes in an intergenic 
region located 50 kb upstream of the Myh2 gene (Figures lA and 
SI) and 4 kb upstream of a lincRNA (2310065F04Rik); we refer to 
this lincRNA as linc-MYH (Figures lA and S2). Sixl binding at 
these MEF3 sites was demonstrated in vivo by ChIP (Chromatin 
Immunoprecipitation) experiments with Sixl antibodies on adult 
fast gastrocnemius plantaris (GP) and tibialis anterior (TA) muscles 
(Figure IB) but not on adult slow Soleus (data not shown), and 
confirmed for five of these sites (sites 1,2,3,4 and 6) by EMSA 
assays (Figure S3A). We asked whether this Myh intergenic region 
could constitute an enhancer element, controlling the spatio- 
temporal expression of Myh genes in this locus. A 2 kb DNA 
fragment of this region, including the six identified MEF3 sites and 
1 kb of DNA fragments upstream of the transcription start site of 
fast-type Myh2, My hi and Myh4 genes, was isolated. The putative 
enhancer was ligated to each Myh promoter using luciferase pGL3 
basic plasmids to generate pGL3-Enhancer-Af);Ai?/7/^ constructs. 
To test the involvement of Six binding in enhancer activation of 
the Myh2, Myhl and Myh4 promoters, we mutated all six MEF3 
sites present in the enhancer, and named these reporters pGL3- 
mutEnhancer-A/jjA^/i/^. Luciferase activity was tested two weeks 
after the electroporation of reporter plasmids in adult TA muscles. 
The luciferase activity of pGL3-Enhancer-ili)'A2/7/^ was seven to 
twelve times higher for either of the promoters, than with pGL3- 
Myh2lll4. Enhancer activity was not observed in plasmids with 



mutated MEF3 sites associated with either of the Myh promoters 
(Figure IC). Enhancer activity was neither observed with the 
promoters of the slow Sin (Figure S3B) or Tnnil genes, or with the 
promoter of the ubiquitous fi-actin gene (Figure S3C). A weak 
enhancer activity was observed with Myh4 promoter in primary 
embryonic fibroblasts, in which Sixl is expressed (Figure S3D and 
data not shown). These data showed that high MYH enhancer 
activity was only reached in vivo and required specific interactions 
with MYH promoter elements. To determine in vivo interactions 
between the enhancer and each Myh gene, we performed 
chromatin conformation capture (3C) assays of adult fast EDL 
(Extensor digitorum longus) myofibers. These experiments re- 
vealed that the enhancer interacts with the promoter of Myh2l 1 / 4 
genes in native chromatin of EDL myonuclei (Figure ID). The 
strongest interactions were observed with the Myhl and Myh4 
promoters, consistent with the expression profile of these two genes 
in EDL muscles. The data demonstrates that the identified 
conserved cis-element acts as an enhancer for the Myh locus and 
that MEF3 sites are essential for its enhancer activity in vivo. 

Loss of S/x7 impairs fast muscle genes and linc-MYH 
expression during post-natal development 

To further characterize the role of Sixl in the control of fast Myh 
gene expression, we bred Sixl^'^^' mice with transgenic mice 
expressing CRE recombinase under the control of the human 
skeletal actin (HSA) promoter and obtained Sixl^"'-^' ;HSA-CRE 
conditional knockout mice (hereafter named cSixl KO) [25,26]. We 
analyzed the expression of fiber type specific genes in the back 
muscles of wild-type control mice and cSixl KO mice at embryonic 
day 18.5 (El 8.5) and at several post-natal stages (two weeks (P2W), 
four weeks (P4W) and eight weeks (P8W)) animals (Figure 2), as 
muscle fiber fast-subtype specialization takes place from the end of 
embryogenesis [9] . Sixl mRNA was not detectable in back muscles 
of cSixl KO mice (Figure 2). The expression of fast-type genes 
{Myh4, Tnnt3, Tnni2, Tnnc2 and Pvalb) increased during postnatal 
development in control mice but that of slow- type genes (Myh7, 
Tnntl, Tnnil, Tnncl and Slri) decreased. The linc-MYH RNA was 
detected after birth in muscle samples and its expression increased 
in line with that of Myh4 (Figure 2). The induction of fast-type 
genes and linc-MYH and the suppression of slow-type genes, were 
impaired in cSixl KO mice. Expression of linc-MYH wsls reduced by 
three to five times in cSixl KO mice during postnatal development 
(Figure 2). These results show that Sixl controls the induction of 
linc-MYH and fast-type genes during postnatal development, and is 
required for the downregulation of slow type genes. 

Sixl deficiency impairs adult muscle fast phenotype 

We next analyzed adult 12 week-old cSixl KO mice to further 
characterize the role of Sixl in adult muscle. Sixl mRNA and 
protein were not detectable in GP enriched with fast-myofibers or 
soleus (SOL) muscle enriched with slow-myofibers (Figure 3A and 
B), and fatigue resistance of TA muscle was 35% higher 
(Figure 30) in the cSixl KO mice. We used immunohistochemistry 
to analyse the composition of MYH7, MYH2 and MYH4 in cSixl 
mutant myofibers. Mutant TA muscles had a higher percentage of 
fibers containing MYH 7 and MYH2, but a lower percentage of 
fibers containing MYH4 (Figures 3D and S4). We found consistent 
results during qPCR analysis oi Myh mRNA i.e., higher levels of 
Myhl and Myh2 mRNA and lower levels of Myh4: mRNA levels 
were observed in the fast TA muscles of cSixl KO (Figure 3E), 
Expression levels of other specific fast and slow-type genes were 
also tested. We found in mutant TA muscles a downregulation of 
fast-type genes {TnntS, Tnni2, Tnnc2 and Pvalb) and a five and to 25 
fold increase in the levels of slow- type genes {Myh7, Tnntl, Tnnil, 
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Figure 1 . Six7 binds directly the enliancer of the fast Myh gene cluster. (A) Schematic representation of the fast Myh gene cluster. (B) qPCR 
values of Chip experiments performed with Sixl antibodies, or IgG on GP and TA chromatin, showing Sixl binding to PI and P2 regions of the fast 
Myh enhancer, to the muscle promoter of AldoA, and to an intergenic region located 86 kb 5' upstream of Myh 3 (NC) (n = 3). (C) Luciferase assays 
from adult TA muscles electroporated with luciferase vectors (indicated) and a TK-Rennila luciferase vector allowing normalization (n = 4). (D) qPCR 
experiments from 3C assays of wild type EDL muscle, showing the direct interactions of Myh2, Myh 1 and Myh4 promoters with the fast Myh enhancer. 
*P<0.05, **P<0.01, ♦**P<0.001. 
doi:1 0.1 371 /journal.pgen.1 004386.g001 



Tnncl and Sin) (Figure 3E). Nevertheless expression of slow Afi'/i7 is 
increased more tlian ten fold at tlie mRNA level in cSixl mutant 
TA myofibers, while by imniunohistochemistry the number of 
MYH 7 positive myofibers is increased less than two fold. This 
showed that there is no major phenotype switch in cSixl mutant 
TA myofibers. This observation could be explained either by the 
higher amount oi Myh7 mRNA accumulating in MYH7 positive 
fibers or by a general increase of Myh 7 mRNA in TA myofibers, 
mRNAs that would not be translated efficiently and leading to the 
absence of increase of MYH 7 positive fibers. The expression of 
linc-MYH expression was lower in the adult TA of cSixl KO mice, 
than in control mice (Figure 3E). These results indicate that the 
Sixl homeoprotein can control the phenotype of fast skeletal 
myofibers in adult animals. 

linc-MYH is expressed exclusively in adult fast-type 
muscles 

We found that linc-MYH is expressed in fast-type skeletal 
muscles (GP, TA and EDL), but not in SOL, brain, kidney, heart 
or fat tissues, an expression pattern which parallels that of the fast- 
fiber Myh4 (Figure 4A). This suggested that linc-MYH is only 
expressed following robust nuclear accumulation of Six 1 , as takes 
place in the nuclei oi MYH4 myofibers [1 1], and that the weaker 
nuclear accumulation of Sixl observed in SOL myonuclei does 
not allow efficient Sixl binding on the MYH enhancer and linc- 
MYH expression. We used luciferase reporter transfection assays 



(as described above) to test the requirement for Six binding on the 
MYH enhancer to activate linc-MYH expression. These transient 
transfection assays, performed in adult TA, show that the MYH 
enhancer activates linc-MYH expression in a Six-dependent 
manner, as measured two weeks after electroporation (Figure 4B). 
lincRMAs can localize in cytoplasm [16] or as a single focus [19] or 
multiple foci [20] in nuclei. To analyze linc-MYH localization in 
skeletal muscle fiber, we performed fluorescent in situ hybridization 
(FISH), using linc-MYH sense and antisense RNA on isolated 
myofibers from fast EDL. Intranuclear localization of linc-MYH 
was observed with the antisense linc-MYH RNA probe, with 
approximately 10 linc-MYH foci per nucleus (n = 10, Figure 4C), 
while the sense linc-MYH RNA probe gave no signal (data not 
shown). We concluded from these experiments that linc-Myh RNA 
accumulates at specific sites in the nucleus of fast myofibers. 

linc-MYH coordinates fiber-type gene expression 

Due to the number of linc-MYH foci observed in fast type nuclei, 
we hypothesized that linc-MYH could act in trans [17] to control 
gene expression in fast myofibers. To test this theory, we used 
electroporation to introduce a shRNA against linc-MYH (shlinc- 
MYH) in TA muscle and analyzed the transfected samples after 
fourteen days. This method yielded the efficient knockdown of linc- 
MYH, with a 90% reduction of its expression (Figure 5A). To 
identify the consequences of linc-MYH knockdown and understand 
its mode of action, RNA samples from shlinc-MYH transfected 
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Figure 2. The expression of fast-type genes and linc-MYH is impaired in cSix7 KO mice during postnatal development. mRNA 
expression level of S/x/, linc-MYH, Myh3, slow-type genes (red) and fast-type genes (blue) In back muscles of cSixl KO mice at El 8.5, P2W, P4W and 
P8W, as determined by qPCR experiments, (n = 3 to 6 for each point). ♦P<0.05, **P<0.01, ***P<0.001. 
doi:1 0.1 371 /journal.pgen.1 004386.g002 



adult TA were analyzed by Afiymetrix microarrays (Figure 5D 
and Table SI), and validated by qPCR experiments (Figure 5A). 
The expression of linc-MYH was significandy lower in the absence 
oi Sixl, but Sixl expression was not afiected by the absence of linc- 
MTH. Knockdown of linc-MYH led to robust gene expression 
modification; this knockdown strongly upregulated the expression 
of numerous slow genes (such as Sin, Tnnil, Tnncl and Tnntl) and 
moderately downregulated the expression of several fast genes, 
(including Mjh4, Tnnt3, Tnni2 and Pvalb) (Figure 5A). We noted 
that the expression of Sin and other slow genes remained far lower 
in /mc-MW/ knockdown TA than in wild type Soleus. This showed 
that linc-MYH knock down in TA is not sufficient to achieve the 
full transcription efficiency of these slow genes as is observed in 
SOL, where calcineurin and several kinases are required to 
achieve their efficient transcription [27]. In addition, and contrary 
to what was observed previously in the muscles of cSixl KO mice, 
slow Mjhl expression level did not change in linc-MYH knock 
down TA. We next tested whether the transcription rate of Sin and 
Tnntl slow genes was modified in linc-MYH knock down TA by 
measuring their pre-mRNA expression level. As can be seen in 
Figure 5B we observed an upregulation of pre-mRNA of these 
slow genes after linc-MYH knock down, which is proportional to 



their respective mRNA accumulation, demonstrating that their 
transcription is increased after linc-MYH knock down in TA, and 
showing that linc-Myh can act in trans to decrease their 
transcription. Since we also observed a moderate down-regulation 
of the expression of fast genes in linc-MYH knock down TA, among 
which Myh4, we next tested whether the activity of the pGL3- 
Enhancer-AfvA^ reporter could be modulated by the absence of 
linc-MYH. Two weeks after cotransfection of pGL3-Enhancer- 
Myh4 and shlinc-MYH in TA, we observed a statistically significant 
decrease of Luciferase activity, demonstrating the requirement of 
linc-MYH for eflTicient transcription of MYH4 (Figure 5C). Alto- 
gether these experiments show that accumulation of linc-Myh in the 
nuclei of fast myofibers participates in the transcriptional silencing 
of slow genes and is required for the fuU activation of fast genes. 

Transcriptomic analysis of cSixl and linc-MYH 

We compared the networks of genes under the control of linc- 
MYH and of Sixl homeoprotein in adult muscles by the 
transcriptomic analysis of cSixl and linc-MYH knockdown 
(Figure 5D and Table S2). We found that the six genes whose 
expression was the most increased in the linc-MYH knockdown 
were also significantly upregulated in cSixl KO muscles (Figure S5). 
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Figure 3. S/x 7 deficiency impairs the adult phenotype of fast muscle. (A) Sixl mRNA expression levels in GP and Sol muscles of three month- 
old control (Ctrl, n = 4) and cSixl KO (n = 3) mice. (B) Western blot analysis of Sixl and ptubulin expression in Sol and GP of Ctrl and cSixl KO mice. (C) 
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Besides slow muscle genes, two genes, Ankrdl and PeglO, were 
more severely upregulated in the linc-MYH knockdown line (10 
and 8 times, respectively) than in cSIXl mutant myofibers (by 2.8 
and 1.5 times, respectively). These non slow-type genes could be 
exclusively repressed by linc-MYH in adult fast myofibers since 
there is a stronger downregulation of linc-MYH accumulation after 
its knockdown than in cSixl mutant myofibers. Transcriptomic 
analysis of adult myofibers deprived of either Sixl or of \inc-Mjh 
identified a strong qualitative and quantitative correlation in the 
expression of specific genes between linc-MYH knockdown and 
cSixl adult mutant myofibers. The expression of slow muscle genes 
was 3 to 1 0 fold higher in linc-MYH knockdown samples, and 5 to 
25 fold higher in C&7KO samples, than in the wild type. We 
further showed that linc-MYH lies downstream of Sixl in the Six 
myogenic pathway and helps to repress slow muscle genes in fast 
myofibers. The downregulation of aU fast-type genes (other than 
Mj>h4:), and the upregulation of slow-type genes, was weaker in the 
linc-MYH knockdown than in the SixlcKO line. Sixl may control 
several inhibitory pathways, including the linc-MYH pathway, to 
prevent slow-type genes expression in adult fast myofibers. During 
fetal development, at a stage where linc-MYH expression is not yet 
activated, Sixl / 4 increases the nuclear accumulation of the slow 
muscle repressors Sox6 and HDAC4 to repress slow muscle gene 
expression [9,28,29]. In accordance with this, the expression of the 
slow genes Myh7, Sin, Tnnil and Tnntl is upregulated in the 



muscle-specific cSox6 mutant [30]. This demonstrates that linc- 
MYH and Sox6, lying both downstream of Sixl, directly 
participate in the downregulation of Sin, Tnnil and Tnntl in fast 
myofibers. However, the repression of slow Mjh7 in fast myofibers 
acts by a Sixl-Sox6 dependent, but linc-MYH independent, 
repression mechanism. In this study, we observed that the levels 
of fast muscle gene expression decreased by 2-3 folds in the linc- 
MYH knockdown, with the highest decrease found for MjM 
expression. The expression of these genes decreased by a factor of 
1.3 to 2.5 in cSixlYJ^, with the highest decrease found for Pvalb 
expression (Figures 3E and 5A). The presence of Six4 and Srx5 
proteins in adult myofibers [11], which have the same DNA 
binding specificity as Sixl, could compensate its absence in cSixl 
KO animals and enable the activation of downstream fast muscle 
targets. In this case, linc-MYH expression could be preferentisiUy 
dependent upon Sixl, rather than on Srx4 or Six5. Altogether, 
these experiments suggest that the accumulation of linc-MYH 
transcripts in the nuclei of fast myofibers facilitates the regulation 
of a network of genes that drive myofiber specialization via the 
same pathway as Sixl and downstream of this transcription factor. 

Forced expression of linc-MYH in adult slow Soleus 

To test the possibility that linc-MYH could modulate the 
expression of specific muscle genes, we forced its expression in 
myogenic C2 cells and in primary myoblasts, where endogenous 
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Figure 4. linc-MYHXs expressed in adult fast-type siteletal muscles and accumulates in their nuclei. (A) Tissue distribution of S/xI, Myh4 
and linc-MYH RNAs. BAT brown adipose tissue; WAT white adipose tissue. (n = 4). (B) Luciferase assays of adult TA electroporated with linc-MYH 
promoter luciferase vectors (Indicated) and a TK-Renllla luciferase vector (allowing normalization). *P<0.05. (C) FISH of Isolated EDL myoflber with a 
linc-MYH antlsense RNA fluorescent-labeled probe (green) and Dapl staining (blue). 
dol:1 0.1 371 /journal.pgen.1 004386.g004 



linc-MYH expression was faintly detectable even in myotubes four 
days after their differentiation (data not shown). Transfection of a 
1 3 kb genomic fragment encompassing the whole linc-MYH gene 
lead to efficient linc-MYH RNA accumulation in myotubes, but 
after this forced expression, we were unable to detect any 
modification in the expression of slow or fast type genes (data 
not shown). These results suggest that specific cofactors of linc- 
MYH required for its appropriate functioning are lacking in 
cultured myotubes in culture, in agreement with the expression of 
linc-Myh only in adult fast type fibers. To circumvent the 
limitations of cultured cells, we turned to in vivo experiments in 
the Soleus in which linc-Myh is weakly expressed. Two weeks after 
linc-MYH gene transfection in the Soleus we observed that linc- 
MYH RNA accumulates up to approximately 80% of its 
expression level in TA (Figure 5E). We observed a selective 
upregulation oi Aijh4 and Pvalb mRNAs which were increased to 
approximately one-third their expression level observed in TA, 
while mRNA for slow genes remained unchanged (Figure 5E). To 
test whether the increase in Mjh4 and Pvalb mRNA was due to an 
increased transcription of their genes, and potentially exclude a 
mechanism implicating mRNA stabilization, we measured pre- 
mRNA accumulation. As can be seen in Figure 5F, the 
transcription of these two fast genes is upregulated in the Soleus 
samples expressing linc-MYH proportionally to their mRNA 
accumulation, demonstrating that linc-Myh can work in trans 
and allow efficient activation of the transcription of specific fast 
genes. Absence of down regulation of the expression of slow genes 
in Soleus myofibers expressing linc-Myh suggests that, like in 
cultured myotubes, linc-Myh RNA needs specific protein-binding 
partners to achieve its function. Such specific protein-binding 
partners may be absent in Soleus as well as in cultured myotubes. 
Nuclear long non coding RNA are known to guide chromatin 



modifiers to specific gene loci, and by recruiting histone modifiers 
or DNA methyltransferase to modulate their transcription rate 
[31]. Potential linc-Myh protein partners expressed differentially in 
fast and slow adult myofibers may explain how linc-Myh efificientiy 
represses the transcription of slow genes and activates the 
transcription of fast genes in fast myofibers, while in slow 
myofibers its forced expression is only able to activate the 
transcription of fast genes. 

Discussion 

The commitment and maintenance of muscle fiber fast sub type 
specialization relies on the specific expression of one of the fast 
Myosin heavy chain gene present at the fast Myh locus, and of 
specific isoforms of sarcomeric genes [1,2,4]. Myosin heavy chains 
are the primary determinant of the efficiency of muscle 
contraction. In this manuscript, we identified a novel mechanism 
for the specialization of the fast-myofiber subtype. We show that 
the long intergenic non-coding RNA linc-MYH and fast MYH 
genes, both of which are essential for myofiber specification, share 
a common enhancer which is regulated by Sixl homeoproteins. 
The linc-MYH specifically accumulates in nuclei of adult fast 
myofibers. Its function, as revealed here by in vivo knockdown and 
transcriptome-wide analysis, is to prevent slow-type muscle gene 
transcription and increase fast-type muscle gene expression in 
fast-type myofibers. We found linc-MYH downregulates the 
transcription of genes associated with slow muscle contractile 
properties like the slow genes Tnn and Sin (a known repressor of 
Sercal/Atp2al protein [32,33] involved in Ca"^ reuptake by the 
sarcoplasmic reticulum). These genes, which belong to the muscle 
contractile machinery and are repressed in adult fast myofiber, are 
positively controlled by Sixl in myogenic C2 cells [34], where linc- 
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Figure 5. Slow-type gene expression is suppressed by linc-MYH. (A) qPCR experiments revealing mRNA expression levels of linc-MYH, Sixl, 
slow-type genes (red) and fast-type genes (blue) in TA muscles expressing a shRNA directed against linc-MYH or LacZ (n = 5). (B) qPCR experiments 
revealing pre-mRNA expression of the Tnntl and Sin slow genes in TA muscles expressing a shRNA directed against linc-MYH or LacZ (n = 5) (C) 
Luciferase activity of pGL3-Enhancer-Myh4 in TA muscles expressing a shRNA directed against linc-MYH or LacZ (n = 4). (D) Microarray analysis of TA 
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doi:1 0.1 371 /journal.pgen.1 004386.g005 
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expression is not detected. This suggests that their expression 
in adult fast myofiber may be restricted by an additional level of 
regulation involving the Sixl-linc-MYH axis. As a result of our 
study we suggest that Sixl controls the acquisition of fast-t>pe 
myofiber mechanical properties by binding to a single enhancer 
region of the fast Myh locus. It promotes the coordinated 
expression of fast Myhs and that of a strong repressor of genes 
controlling slow contractile properties. The modulation of Six 
activity (depending on fiber-t\pe) facilitates changes in the 
expression levels of the fast genes Myh and Tnn; these changes 
are required for the formation of efficient sarcomcric units and the 
appropriate Ca"^ cycling and excitation/ contraction/relaxation 
coupling [1^]. The Myh enhancer element therefore connects 
distinct regulatory hubs to achieve ultimate muscle fiber special- 
ization. In this context, linc-MYH functions as an end-of-the-chain 
control element, conveying information on the state of fast Myh 
enhancer activity to repress slow-type specific genes and coordi- 
nates a finer level of regulation. This genomic organization at the 
fast Myh locus is reminiscent of the slow Myh7 locus where two 
microRNA miR-208b and miR-499 involved in fast myofiber 
program repression are co-regulated with Myh? [35]. The precise 
molecular interactions between linc-Myh and higher order chro- 
matin modifying" complexes remains to be identified, to explain 
how linc-Myh coordinates the activation of target genes at specific 
sites in the nucleus, and the repression of others. 

Materials and Methods 

Mice, ethics statement 

Animals were bred and handled as recommended by European 
Community guidelines. Experiments were performed in accor- 
dance with the guidelines of the French Veterinary Department. 
t&iKO mice were obtained by breeding the Sixl-LoxP mice [26] 
and transgenic mice expressing a CRE recombinase under the 
control of the human skeletal actin promoter (HSA) [25] . 

Chip experiments 

GP and TA muscles of 2 months old mice were minced with 
scissors just after sampling and fixed in 1 % formaldehyde for 
10 minutes. Formaldehyde was quenched by addition of 0.125 M 
glycine, and muscles were washed twice in PBS. The muscles were 
incubated on ice in lysis buffer (10 mM Tris-HCI pH 7.9, 85 mM 
KCl, 0.5% NP40, protease inhibitors (cOmplete, Roche)) for 
10 minutes and homogenized with a mortar and, subsequendy 
with a Dounce homogenizer. The nuclei were obtained by 
centrifugation, incubated in SDS lysis buffer (50 mM Tris-HCI 
pH 8, 10 mM EDTA, 1% SDS, protease inhibitors) for 10 min- 
utes, and sonicated in a Bioruptor apparatus (Diagenode). The 
debris were removed by centrifugation. The sonicated DNA was 
incubated with 1 [ig of Sixl antibodies (HPAOO 189.3, Sigma) under 
agitation at +4°C overnight. 20 |i,l of Dynabeads protein G 
(Invitrogen) were added to the samples and incubated under 
rotation at +4°C for 1 hour. The beads were washed with low-salt 
buffer (2 mM EDTA, 20 mM Tris-HCI pH 8, 150 mM NaCl, 1% 
TritonX-100, 0.1% SDS), high salt buffer (2 mM EDTA, 20 mM 
Tris-HCI pH 8, 0.5M NaCl, 1% TritonX-100, 0.1% SDS), LiCl 
buffer (1 mM EDTA, 10 mM Tris-HCI pH 8, 0.25M LiCl, 1% 
NP40, 1% deoxycholate) and TE buffer (1 mM EDTA, 10 mM 
Tris-HCI pH 8).' The DNA was eluted with elution buffer (1% 
SDS, 0.1 M NaHCOs) containing 0.1 mg/ml proteinase K 
(Invitrogen) at 62°C for 2 hours, and, proteinase K was inactivated 
by incubation at 95°C for 10 minutes. The DNA was finally 
purified with MinElute PCR purification kit (Qiagen). The 
amount of specific amplified DNA is normalized by beta-Actin 



promoter amplification. The sequences of the oligonucleotides 
used in this study are as foUows. Enh IF, 5'-ATC TCC ACC TCC 
CTC CAA CT; Enh IR, 5'-ACC CCC TAG CTT TGA CAG 

GT; Enh 2F, 5'-AAT CTG ACG ACA GGG TGA GC; Enh 2R, 
5'-GGT CGC CTG ACC TGA TAG AG; AldoaF, 5'-CTC TCA 
AGG CAA ACC AAA GC; AldoaR, 5'-CCA GTG TCC CAG 
ACC TTC TC; ActbF, 5'-TGT TAC CAA CTG GGA CGA CA; 
ActbR, 5'-ACC TGG GTC ATC TTT TCA CG, NCF, 5'-ATC 
CTG CCC CAC TGT GTT AG; NCR, 5'-GCC AGC AAT 
TTG GTT TGA AT. 

3C experiments 

3C experiments were performed as described [36] with few 
modifications. Single myofibers were obtained from adult EDL 
muscles as previously described [26], cross-linked in 2% formal- 
dehyde, 10 mM Tris-HCI pH 7.9, 85 mM KCl, 0.5% NP40 for 
10 min at room temperature. Crosslinking reaction was quenched 
by 1 M glycine. Cross-linked myofibers were lysed for 10 minutes 
with lysis buffer (10 mM Tris-HCI pH 7.9, 85 mM KCl, 0.5% 
NP40) on ice, and the nuclei were harvested. Nuclei were 
resuspended in appropriate restriction enzyme buffer, 0.3% SDS 
and incubated for 1 hour at 37°C with shaking. Triton X-100 was 
added to 2%, and samples were incubated for 1 hour at 37°C. 
Samples were digested with Hind III overnight at 37°C. DNA 
ligation was performed for 4 hours at 16°C and for 30 minutes at 
room temperature". Cross-links were rcv(;rsed, and DNA was then 
purified by phenol extraction and ethanol precipitation. To correct 
for the PCR amplification efficiency of differ(-nt primer sets, a 
BAC clone containing the mouse Myh locus (RP23-61C14) was 
digested, ligated and used as control templates. Quantification of 
the data was performed by quantitative real-time PCR using the 
Lightcycler 480 probe master (Roche Diagnostic). The sequences 
of the oligonucleotides used in this study are given in Table S3. 

RNA preparation 

TA, back, soleus and GP muscles were collected from cSixl KO 
and control mice. Total RNAs were extracted by Trizol Reagent 
(Invitrogen) according to manufacturer's instruction. 

cDNA synthesis and qPCR 

RNAs were treated with DNase I (Turbo DNA-free, Invitrogen) 
and were reverse transcribed with Superscript III kit (Invitrogen) 
according to manufacture's instruction. Reverse transcription was 
performed with 1 |ig of total RNA. Quantitative real time PCR 
(Light Cycler 480, Roche) was performed using Light Cycler 480 
SYBR Green I Master Kit (Roche) according to the manufactur- 
er's protocols. PCR was performed for 40 cycles of 95°C for 
15 seconds, 60°C for 15 seconds, and 72°C for 15 seconds. Genes 
expression level was normalized by the expression level of the 
housekeeping gene Actb. The sequences of the oligonucleotides 
used in this study are given in Table S4. Pre-mRNA qPCR 
experiments to measure RNA transcription rate were performed in 
the same conditions. Reverse oligonucleotides were complemen- 
tary to intronic sequences, while forward oUgonucleotides were 
complementary to exonic sequences. Samples without reverse 
transcription were used as controls, and signal due to contam- 
inating DNA was subtracted to the values obtained with cDNA. 
We noticed that genomic DNA contamination was very low (less 
than hundred fold level of qPCR value observed with cDNA). 

Muscle contraction test 

Skeletal muscle function was evaluated by measuring in situ 
muscle contraction, as described previously [37]. 12 week-old male 
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mice were anesthetized (intraperitoneal injection of pentobarbital 
sodium, 50 mg/kg). Body temperature was maintained at 37°C 
using radiant heat. The distal tendon of the TA muscle was 
attached to an isometric transducer (Harvard Bioscience) using a 
silk ligature. The sciatic nerves were proximally crushed and 
distally stimulated by a bipolar silver electrode using supramax- 
imal square wave pulses of 0. 1 ms duration. Responses to tetanic 
stimulation (pulse frequency 50-143 Hz) were successively record- 
ed. Maximal forces were determined at optimal length (length at 
which maximal force was obtained during the tetanus). Fatigue 
resistance was then determined after a 5-minutes rest period. The 
muscle was continuously stimulated at 50 Hz for 2 minutes (sub- 
maximal continuous tetanus), and the duration corresponding to a 
20% decrease in force was recorded. 

RNA-FISH 

Fluorescent-labeled antisense linc-MTH probes were synthesized 
according to manufacturer's instruction (FISH Tag RNA kit, 

Invitrogen). FISH experiments were performed on isolated EDL 
myofibres and images acquired on a Leica SP2 confocal 
microscope. 

Generation of shRNA against mouse linc-MYH 

Five distinct shRNAs targeting mouse linc-MYH were designed, 
called MincMTH, and inserted into the psiSTRIKE hMGFP 
system (Promega). The efficiency of each shRNA was established 
by determination of linc-MTH transcript levels in TA muscles 
transfected by each shlincMTH. The shRNA against 5'-TTC TGC 
TCA CCA CCT AC A ATT-3' sequence was selected for the 
knockdown experiment. For knock down experiments using 
shlincMTH, a plasmid coding for shLac^ was electroporated in 
the contra-lateral TA as a negative control. 

Eiectroporatlon 

In vivo transfections were also carried out on 10-weeks old 
C57B16 mice. For each experimental conditions, three to five 
Tibialis anterior (TA) or Soleus (Sol) muscles belonging to different 
mice were used. Under isoflurane anesthesia, legs were shaved and 
muscles were pre-treated by injection of a sterile 0.9% NaCl 
solution containing 0.4 U of bovine hyaluronidase/|J,l two hours 
before plasmid injection. Ten |.ig of shRNA-expressing vector were 
introduced into TA muscles of 8 week-old mice by eiectroporatlon 
as previously described [11]. Two weeks following electroporation, 
TA myofibers expressing GFP were dissected under a Nikon 
SMZ1500 stereo microscope and frozen in liquid nitrogen before 
processing for Luciferase assays or RNA purification. 

Immunohlstochemistry 

TA, soleus and gastrocnemius muscles were embedded in 
cryomatrix and quickly frozen in isopentane cooled with liquid 
nitrogen. Cryostat sections (10 |j.m) were fixed in 4% PFA, washed 
in PBS, permeabUized with 0.1% Triton X- 1 00 and left for 1 hour 
in blocking solution (1 X PBS, 1.5% goat serum, 0.1% Triton X- 
100). Rabbit poly-clonal antibodies directed against Laminin 
(Z0097, Dako) (1/100 dilution), and monoclonal antibodies 
against MTH7 (NOQ,7.5.4D, Sigma) (1/1000 dilution), MTH2 
(SC-71, Developmental Studies Hybridoma Bank) (1/20 dilution) 
and against MTH'i (BF-F3, Developmental Studies Hybridoma 
Bank) (1/20 dilution) were applied overnight at 4°C to the treated 
sections. The next day, after three washes with 1 x PBS containing 
0.05% Tween-20, cryosections were incubated for 1 h with 
appropriate ffuorescent secondary antibodies (Alexa Fluor 488 
goat anti-rabbit IgG 1/1000 dilution, Alexa Fluor 594 goat 



anti-mouse IgG 1/1000 dilution, Invitrogen). After three washes 
with 1 X PBS containing 0.05% Tween 20, samples were mounted 
in Vectashield mounting medium. 

Microarray 

After validation of RNA quality with the Bioanalyzer 2100 
(using Agilent RNA6000 nano chip kit), 50 ng of total RNA were 
reverse transcribed following the Ovation PicoSL WTA System 
(Nugen). Brieffy, the resulting double-strand cDNA was used for 
amplification based on SPIA technology. After purification 
according to Nugen protocol, 5 (Xg of single strand DNA was 
used for generation of Sens Target DNA using Ovation Exon 
Module kit (Nugen). 2.5 |J,g of Sens Target DNA were fragmented 
and labelled with biotin using Encore Biotin Module kit (Nugen). 
After control of fragmentation using Bioanalyzer 2100, the cDNA 
was then hybridized to GeneChip Mouse Gene 1.0 ST (Affyme- 
trix) at 45°C for 17 hours. After overnight hybridization, the 
ChlPs were washed using the fluidic station FS450 following 
specific protocols (Affymetrix) and scanned using the GCS3000 
7G. The scanned images were then analyzed with Expression 
Console software (Affymetrix) to obtain raw data (eel files) and 
metrics for Quality Controls. The analysis of some of these metrics 
and the study of the distribution of raw data show no outlier 
experiment. RMA normalization was performed using R and 
normahzed data was subjected to statistical tests. 

EMSA 

EMSA was carried out with Sixl full-length mouse cDNA 
cloned into the pCR3 vector (Clontech) as previously described 
[38]. Recombinant mouse Sixl protein was obtained with a T7 
transcription/translation kit (Promega). The oligonucleotide con- 
taining double-stranded myogenin MEF3 site was incubated with 
recombinant proteins. Competition experiments were performed 
in the presence of a ten-fold and hundred-fold molar excess of 
unlabeled identified Myh enhancer MEF3 sites (Enhl to Enh6) or 
mutated A4yh MEF3 sites (mtEnhl to mtEnhS), or Myogenin 
promoter NFI or MEF3 sites. The sequences of the oligonucle- 
otides used are as follows, the MEF3 consensus sequence is 
underlined; EnhlF 5'-CTC TTG GGT AAC TGG AGC CCC 
TC-3'. EnhlR 5'-GAG GGG C TC CAG TTA C CC AAG AG- 
3'. Enh2R 5'-GGT TGA C TT AGA TTT C CT TAT GA-3'. 
Enh2F 5'-TCA TAA GGA AAT CTA ACT CAA CC-3'. Enh3F 
5'-TGT AAG AGA AAC TGA AAT AAA AT-3'. EnhSR 5'- 
ATT TTA TT T CAG TTT C TC TTA CA-3'. Enh4F 5'-GGG 
GTA AGA AAT CTG ACG ACA GG-3'. Enh4R 5'-CCT GTC 
G TC AGA TTT C TT ACC CC-3'. EnhSF 5'-CTA TCA GGT 
CAGGCGACC TCA GT-3'. Enh5R 5'-ACT GAG GTC GCC 
TGA CCT GAT AG-3'. Enh6F 5'-CGT CAA GGA AAC CTT 
ATT CCA TC-3'. Enh6R 5'-GAT GGA A TA AGG TTT C CT 
TGA CG-3'. MyogF 5'-TGG GGG GGC TCA GGT TTC TGT 
GGC GT-3'. MyogR 5'-ACG CCA CAG AAA CCT GAG CCC 
CCC CA-3'. NFIF 5'-TAT CTC TGG GTT CAT GCC AGC 
AGG G-3'. NFIR 5'-CCC TGC TGG CAT GAA CCC AGA 
GAT A-3'. mtEnhlF 5'-CTC TTG GGT AGG ATC GGC CCC 
TC-3'. mtEnhlR 5'-GAG GGG CGG ATC CTA CCC AAG 
AG-3'. mtEnh2F 5'-GGT TGA CGA ATT CTT GCT TAT 
GA-3'. mtEnh2R 5'-TCA TAA GCA AGA ATT CGT CAA CC- 
3'. mtEnh3F 5'-TGT AAG ACC AAC TGA AAT AAA AT-3'. 
mtEnh3R 5'-ATT TTA TTT CAG TTG GTC TTA CA-3'. 
mtEnh4F 5' -GGG GTA AGA AGG ATC CCG ACA GG-3'. 
mtEnh4R 5'-CCT GTC GGG ATC CTT CTT ACC CC-3'. 
mtEnhSF 5'-CTA TCA GGT CGG ATC CCC TCA GT-3'. 
mtEnh5R 5'-ACT GAG GGG ATC CGA CCT GAT AG-3'. 
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mtEnhSF 5'-CGT CAA GGA AGO ATC CTT CCA TC-3'. 
mtEnh6R 5'-GAT GGA AGG ATC CTT CCT TGA CG-3'. 

Western blot 

Western blots were performed with protein extracts of GP and 
soleus muscles from c.Su7KO mice and control mice as previ- 
ously described [9]. 1:1000 dilutions of anti-Six 1 antibodies 
(HPA001893, Sigma) or anti-p- tubulin antibodies (2128, CeU 
Signaling) were used. 

Statistical analysis 

All graphs represent mean values ± SEM. Significant differ- 
ences between mean values were evaluated using two-tailed, 
unpaired Student's t test (when two groups were analyzed) or one- 
way ANOVA followed by Student Newman-Keuls test (for three 
or more groups). 

Plasmids construction 

For the construction of the pGL?>-Mjh2/l /4, c57bl6N mouse 
DNA was first used as a template to clone 1 . 1 kbp promoters of 
Myh2/l/4 with forward KpnI/SacI 5'-TTC AGA AAC TGC 
ATC ACT TAA A-3' and reverse Mlul, 5'-GCA GCT CGG 
GCA GTG GCC AGT GT-3', forward KpnI/SacI 5'- CAT 
ATC TGC ATC TCT AGA TAC C-3' and reverse Mlul, 5'- 
GGC AGC AGC AGC CAG GAT GTG T-3', forward Kpnl/ 
Sad 5'- ACC GCT AGC CTT GAG CCT TTG-3' and reverse 
Mlul, 5'- ATA GCG AGA GCC CTT TGT TCT C-3', 
respectively. Myh2/l/4 promoter fragments were subsequendy 
inserted into an KpnI-MluI digested pGL3 basic plasmid. For the 
construction of the pGL3-Enhancer-A/))A2/i/^, mouse DNA was 
first used as a template to clone the enhancer with forward Kpnl 
5'- GCG TTT CTA ATT CGG CTT GAA C-3' and reverse 
Sad, 5'- CAT TTC CTT CCT CTA AAG GCT CTT TATT C- 
3'. This enhancer fragment was subsequently inserted into Kpnl- 
SacI digested pGL'i-Myh2 / 1 / 4 plasmids. For the construction of 
the pGL3-mutEnhancer-A^A2/i/^, the six MEF3 sites of the 
enhancer were mutated as follows; MEF3 1: 5'-GTAACTGGA to 
5'-GTAGGATCC; MEF3 2: 5'-TTAGATTTC to 5'-GAAT- 
TCTTG; MEF3 3: 5'-GAAACTGAA to 5'-CCAACTGAA; 
MEF3 4: 5'-GAAATCTGA to 5'-GAAGGATCC; MEF3 5: 5'- 
GTCAGGCGA to 5'-GTCGGATCC; MEF3 6: 5'-GAAA- 
CCTTA to 5'-GAAGGATCC. All plasmids sequence was 
confirmed by sequencing. 

For the construction of the pSF-pA-CMVe-/mc-iV/17/, genomic 
DNA fragment containing linc-MYH was obtained from digestion 
of a BAG clone containing (RP23-61C14) by BsaWI and Avrll. 
The 13.3 kbp DNA fragment was subsequently inserted into a 
Spel-Xmal digested pSF-pA-CMVe plasmid. For linc-Myh gain of 
function experiments, the empty pSF-pA-CMVe plasmid was 
electroporated in the contra-lateral Soleus as a negative control. 
For the construction of the -pGLS-Adb, pGL3-6'/n and pCjh'i-Tnnil , 
mouse DNA was first used as a template to clone the promoters of 
Actb, Sin, and Tmil with forward 5'- TTT CTC TAT GGA TAG 
GTA CCT TTG AGC TCC TGA CCC CGT GTG TAG CTC 
T-3' and reverse 5'- GAG CCC GGG CTA GCA CGC GTA 
AGG AGC TGC AAA GAA GCT G-3', forward 5'- TTT CTC 
TAT CGA TAG GTA CCT TTG AGC TCT ACC GAC TAT 
CAT GCC CAC A-3' and reverse Mlul, 5'- GAG CCC GGG 
CTA GCA CGC GTC AGG CTA CCA AGG ACC TCA G-3', 
forward 5'- TTT CTC TAT CGA TAG GTA CCT TTG AGC 
TCC TGG GAT TTG AAC CCA TGA C-3' and reverse 5'- 
GAG CCC GGG CTA GCA CGC GTC CTC ACC ACA GAC 
TGC AGA G-3', respectively. Actb, Sin, and Tnnil promoter 
fragments were subsequendy inserted into a KpnI-MluI site of 



pGL3 basic plasmid by GeneArt kit (Life Technologies). For the 
construction of the pGL3-Enhancer-^riA, pGL3-Enhancer-iS'/n and 
pGL3-Enhancer-T«nz7, the enhancer fragment was subsequendy 
inserted into an KpnI-SacI site of pGL,3-Actb, pGL,3-Sln and 
pGL3-T?!«27plasmids. AU plasmids sequence were confirmed by 
sequencing. 

Luciferase assays 

Two |J.g of Luciferase-expressing vector and one hundred ng of 
pRL-TK vector (Promega) were introduced into TA muscles of 8 
week-old mice by electroporation as previously described [II]. 
Two weeks following electroporation, the TA muscles were 
dissected and frozen in liquid nitrogen before processing. The 
TA muscles were homogenized in Passive Lysis Buffer (Dual- 
Luciferase Reporter Assay System, Promega) and rotated for 
15 minutes. The homogenate were centrifuged to remove debris, 
and the supernatant was used for measurement of Luciferase 
activity according to manufacture's instruction (Dual-Luciferase 
Reporter Assay System, Promega). 

Computational analysis 

In order to computationally identify MEF3 binding sites, we 
buUt a PWM (position-specific weight matrix) for MEF3, starting 
from a list of 15 binding sites (see Table S5) that were previously 
tested by Electrophoretic Mobility Shift Assay on the basis of their 
proximity to the MEF3 Myogenin consensus GAAACCTGA [10]. 
The PWM (shown in Figure S6 and Table S6) was generated by 
running the de novo motif finder Imogene [39] on small DNA 
fragments (see Table S5) that contained these binding sites. 
Imogene used phylogeny to enrich mouse set of DNA fragments 
with orthologs in 1 1 other mammalian sequenced genomes and to 
produce a refined PWM. The information content of the PWM 
(the genmot Sg parameter of Imogene) was set to 8.7 bits. Binding 
sites were predicted using a prediction threshold (the scangen Ss 
parameter of Imogene) of 9 bits and requiring conservation, as 
explained in [39]. 

Supporting Information 

Figure SI Sequences of PI and P2 boxes of the Alyh enhancer. 
Sequences of PI and P2 boxes of the Myh enhancer in mouse, rat, 
human, bovine and equinides species, and showing the sequence 
conservation of the six MEF3 sites and E boxes. Coordinates are in 
mm9 assembly. 
(TIF) 

Figure S2 Predicted linc-MYH structure. Predicted minimum 
free energy (MFE) structure of the 1050 nt long linc-MYH, as 
determined by RNAfold [40]. The color encodes base-pair 
probabilities. 
(TIF) 

Figure S3 (A). Competitive Electromobility shift assays. Com- 
petitive Electromobility shift assays performed with recombinant 
Sixl proteins and labeled Myogenin MEF3 oligonucleotide and 10 
or 100 fold molar excess of unlabeUed oligonucleotides containing 
Myogenin MEF3 or NFl site, with MYH MEF3 sites (1, 2, 3, 4, 5, 
6) or with mutated MYH MEF3 sites whose sequence is presented 
on Figure SI and in the Materials and Methods section. (B). qPCR 
experiments measuring the relative levels of Sin mRNA in adult wt 
TA and Sol. (C). Luciferase assays of adult TA transfected with Sin, 
Tnnil and Actb promoters with or without the Myh enhancer. TA 
sampling was performed two weeks after the electroporation, 
(n = 4). (D). Primary embryonic fibroblasts were transfected with 
Luciferase plasmids under the control of Myh4 promoter or Myh4 
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promoter linked with the Myh enhancer. After transfection, 
fibroblasts were cultured two days before sampling, (n = 3). *P< 
0.05, ***P<0.001. 
(TIF) 

Figure S4 Immunostaining of MYH proteins. Immunostaining 
of MYH7 (red), MYH2 (red), MYH4 (red) and laminin (green) in 
TA of 12 weeks old control and cSixlKO male mice. 
(TIF) 

Figure S5 Comparison between cSixl mice and sklinc-MTH 

treated mice. (-^4) Venn diagram showing the overlap between 
genes that are up-regulated more than five fold in ihlinc-MYH and 
two fold in cSixl muscles, {p— 10~" as given by a hypergeometric 
test). {B) Scatter plot of mRNA expression fold change (log2) as 
determined by microarray analysis of sh/mc-AfKF/ and cSbcl. Genes 
that are up-regulated more than five fold in shlinc-MTH and two 
fold in cSixl muscles are indicated in red. (C) mRNA expression 
levc-ls in shlinc-AITH knock-down and cSixl TA muscles, as 
measured by qPCR experiments. Ctrl (means ± SEM; n = 4), 
cSixl (means ± SEM; n = 4), shiac^ (means ± SEM; n = 5), shlinc- 
MTH (means ± SEM; n = 5). (Z)) mRNA expression level of 
Ankrdl, Z'^hf2 and PeglO in back muscles of cSbcl KO mice at 
El 8.5, P2W, P4W and P8W, as determined by qPCR experi- 
ments, (means ± SEM; n = 3 to 6 for each point). *P<0.05, **P< 
0.01, ***P<0.001. 
(TIF) 

Figure S6 MEF3 PWM used in diis study. 
(TIF) 

Table SI Microarrays analysis of TA muscles electroporated by 

shRNA. 

PCLS) 



Table S2 Microarrays analysis of GP muscles of cSixl. 

(XLS) 

Table S3 Sequence of the oligonucleotides used for 3C. 
(DOCX) 

Table S4 Sequence of the oligonucleotides used for qPCR. 

(DOCX) 

Table S5 Coordinates (mm9) of the fragments used to learn the 

MEF3 PWM. 

(DOCX) 

Table S6 MEF3 frequency matrix. 
(DOCX) 
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